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Abstract 

The detection of a soft thermal X-ray component in the spectrum of a bright knot in the halo of the plerion G21.5- 
0.9 is reported. Using a collisional ionization equilibrium model for an hot optically thin plasma, a temperature 
kT = 0.12 — 0.24 kev, a mass of 0.3-1.0 M© and a density of 1.6-6 cm -3 is derived. The spectral analysis suggests 
a possible overabundance of Silicon with respect to the solar value in the knot; if this will be confirmed this object 
may be a clump of shocked ejecta. 
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1. Introduction and data analysis 

The plerion G21.5- 0.9 has been the subject 
of various studies in radio ( Bock et al.l200llandref- 
erences th erein), infrared GaUant^md^iffs 1998) 
and X -ray ijSafi-Harb et alJ2CjolIlL£ri^al^rnbara and 
Hqq2). However its weak X-ray hal o was discovered 
only recently ( Slane et all l20f)(i I War wick et al.1 
l200ll ) thanks to the larger sensitivity of Chan- 
dra llWeisskopf et all l2000|) and XMM-Newton 
Hansen et al.l l200l|) with respect to previous X- 



namely the North Spur ((Warwick et al.ll200lh . 

For XMM-Newton, the data consist of the 
observation of G21.5-0.9 with the source posi- 
tioned on the axis and off the axis of the EPIC 
cameras (observation ID 0122700101, 0122700201, 



ray astronomy satellites. The weakness of the halo 
has prevented a detailed study of its features up 
to now. Luckily, G21.5-0.9 is a calibration source 
for both Chandra and XMM-Newton and a large 
number of calibration datasets have been accu- 
mulated and is still being accumulated to this 
purpose. In this paper, this large collection of 
public data is used to investigate the nature of the 
brightest spot in the X-ray halo of this plerion, 



McrQilMPOl, 0122700401, 0122700501), which have 
a total summed live time of 142 ks for PN and 146 
ks for the MOS cameras. The observations were 
taken on April 7-17, 2000. The original event files 
were screened to eliminate the contribution of soft 
proton flares. After the screening, the summed 
exposure time were 79 ks for PN and 116 ks for 
MOS. All the analysis of XMM data was per- 
formed with the software SAS v6.0. The XMM- 
Newton spectra of the North Spur, collected in 
the single observations using a circular region of 
20 arcsec radius centered at location 18' l 33 m 32.7 s 
RA and -10 d 32 m 49 s DEC (J2000), were summed 
and their response averaged. The spectra of the 
two MOS cameras were also summed and their 
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effective areas and response matrices were aver- 
aged. The fit was performed separately for the PN 
and MOS, and then, after checking that they give 
consistent results, they were fitted simultaneously. 
The spectra have been rebinned to a minimum of 
20 counts per bin and the 0.5-9 keV energy range 
was selected for the spectral fitting. The back- 
ground used for the spectral analysis was collected 
m an annulus centered on 18 /l 33 m 33.8 s RA and 
-10 d 34 ro 07 s DEC J2000 between 85" and 120" 
from the center (this region is inside the X-ray halo 
at the same average off-axis angle as the North 
Spur). This method subtracts the contribution of 
the underlying continuum due to diffuse emission 
of the X-ray halo from the source spectrum, and 
concentrate the study on the intrinsic spectral 
properties of the North Spur. 

For Chandra, 21 observations from the 
Chandra archive for which the source was located 
on the ACIS S3 chip were selected, in order to 
maximize the instrument response in the soft X- 
ray band, and at off-axis angle < 5'. The sequence 
number of the selected observation are 0159, 1433, 
1554, 1717, 1769, 1771, 1839, 2873, 3693, 4353, 
5166,1230, 1553, 1716, 1718, 1770, 1838, 1840, 
3474,37 00, 4354. The da t asets include the dataset 
used bv ISafi-Harb et all l|200lf) . the latest report 
on X-ray data of the North Spur. However, they 
used only 6 observation for a total of 65 ks, whereas 
we also include all the more recent observation 
that have been done, for a total of 196.5 ks. This 
is the reason why they did not detect the thermal 
emission reported in this work. The observations 
were reprocessed using CIAO 3.1 and CALDB 
2.28 in order to apply the most recent calibra- 
tions to all the data, including time dependent 
gain correction. Afterwards, the data were filtered 
to retain grade 0,2,3,4,6 and the status column 
equal zero. All the filtered datasets were merged 
together using the CIAO merge_ALL task. The 
spectrum of the North Spur was extracted using 
a circular region centered at the same location as 
before, but with a smaller radius (13 arcsec) to 
exploit the high resolution of the Chandra mir- 
rors. The spectrum was rebinned like the EPIC 
spectrum and the band 0.5-8 keV was used. The 
response matrix and effective area were generated 
using the CIAO tasks for the analysis of extended 



sources. The background was selected in the same 
way as the XMM-Newton background. All the fits 
were performed using XSPEC vll.3. 



2. Results of spectral fitting 

The spectra were first fitted with a power 
la w model modified by the interstel lar absorption 
of lMorrison and McCammonl ( 1983h and a reason- 
able fit was found. All the best-fit parameters and 
their uncertainties are reported in Table 1. How- 
ever, given the non uniform distribution of resid- 
uals in the soft band and the possibility that the 
North Spur may show signs of thermal emission 
due to hot plasma, as part of a shell or ejecta asso- 
ciated to G21.5-0.9, the spectra were fitted with a 
combination of a power law model an d the MEKAL 
optically-thin thermal pl asma model l|Mewe et alJ 
ll985ilLiedahl et al.ll995h . The results, also showed 
in Table 1, indicate that the two components fit 
is better then the power law only fit. The decre- 
ment of x 2 is significant according to the F-test 
(probability ~ 10~ 4 ), and the distribution of resid- 
uals is more homogeneous. The discrepancy be- 
tween XMM-Newton and Chandra derived fluxes 
are due to different extraction region area and sug- 
gests that the source is extended (especially its 
non-thermal emission) . 

The XMM-Newton and Chandra spectra 
are shown in Fig. 1 and 2, along with the best-fit 
two component model. Emission lines of Mg XI 
at 1.34 keV and Si XIII at 1.86 keV may be iden- 
tified in Fig. 2. Other expected emission lines at 
different energies are not evident due to limited 
statistics of these highly absorbed spectra and 
moderate CCD spectral resolution. 

Fig. 3 shows the x 2 contour levels in the 
temperature-emission measure plane correspond- 
ing to confidence levels of 68%, 90% and 99% on pa- 
rameter uncertainties. Fig. 4 shows the same con- 
tours levels of Fig. 3, but for the Mg and Si abun- 
dances when these parameters are left free to vary. 
The x 2 obtained in this case is better than the case 
of abundances fixed to solar values (F-test proba- 
bility of 2 x 10~ 3 ), but the solar values for both 
elemental abundances are still within the 99% con- 
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Table 1 

Summary of spectral fitting results to the X/ray spectrum of the bright knot "North Spur. 



Model 


Nh 7 


Norm 


kT 


flux a 


X 2 /dof 




cm -2 photon cm 


" 2 s" 1 keV 


- 1 keV 


erg cm -2 s~ 


i 




X 10 22 index 


Xl0~ 4 




xlO" 15 






XMM-Newton PN and MOS 






Power-law 


1.45 + 0.13 1.91 + 0.12 


13.9 + 2.5 






442/407 


P.L.+MEKAL 


2.63 + 0.30 2.30 + 0.15 


29.3 + 6.5 


0.15 ± 0.03 


12.0 + 4.5 


423/405 


P.L.+VMEKAL b 


2.63 + 0.39 2.26 + 0.13 


28.5 + 7.6 


0.17 + 0.03 


11.4 + 5.6 


410/403 




Chandra ACIS S3 








Power-law 


1.15 + 0.09 1.72 + 0.10 


9.3 ± 1.2 






237/187 


P.L.+MEKAL 


1.78 + 0.32 2.04 + 0.16 


15.4 + 3.7 


0.17 + 0.04 


5.3 + 1.3 


217/185 


P.L.+VMEKAL c 1.81 ± 0.24 2.01 ± 0.14 


14.8 + 3.0 


0.15 ± 0.03 


7.9 + 2.3 


206/183 


a Absorbed flux due to the thermal component only, in 


the 0.5-2.0 keV. 





b Only Mg and Si abundances were fitted (best-fit values and 90% uncertainties are 1.0(0.1- 
2.3) and 4(1.3-8), respectively). Flux of Mg XI line at 1.34 keV is 3.41"^ x 10~ 6 ph cm" 2 
s~\ whereas flux of Si XIII line at 1.86 keV is 8.7tfl x 10" 6 ph cm" 2 s" 1 . 
c Only Mg and Si abundances were fitted (90% confidence range is > 0.6 for Mg and > 5 for 
Si). Flux of Mg XI line at 1.34 keV is 2.8tl' \ X 10" 6 ph cm" 2 s" 1 , whereas flux of Si XIII 
line at 1.86 keV is 9.5±l'l x 10" 6 ph cm" 2 s" 1 . 



fidcilCC level Chandra ACIS S3 spectrum of the North Spur 

Best-fit VMEKAL+POWLAW model 




with a thermal+non thermal model. We also show the Energy (keV) 

contribution of thermal and non-thermal components to 

the total best-fit spectra (for MOS only). Fig. 2. Same as Fig. 1 but for Chandra ACIS S3. 
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Fig. 3. 68%, 90% and 99% confidence level contours for the 
plasma temperature and emission measure of the thermal 
component used in the fit of the North Spur XMM-Newton 
EPIC (thin lines) and Chandra ACIS-S (thick lines) spec- 
tra. 
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Fig. 4. Same as Fig. 3 but for Mg and Si abundances 
relative to solar. 

3. Discussion 

The best-fit Si abundances obtained from 
the fit is around 4 and its 90% confidence level 
range is between about 2 and 10 for the XMM- 
Newton EPIC and above 5 for Chandra ACIS-S. 
Even though the solar value cannot be excluded 
at a 99% confidence level for XMM-Newton, the 
fact that both Chandra and XMM-Newton inde- 
pendently give a best-fit Si abundance above solar 
can be considered a hint that this element is some- 
what overabundant with respect to solar. On the 



other hand, Mg abundance is found to be compat- 
ible with the solar value. 

The best-fit Si abundance is not compati- 
ble with the abundances found in shocks of young 
SNR inte racting with the circumstellar material 
(see, e.g.. IPark et alJl2004h . Instead, it may indi- 
cate that the North Spur is an ejecta knot. The 
fact that the Mg abundance ranges between 0.3 
and 2 with a best-fit of ~ 1 may imply that we 
are seeing at a Si-rich knot like the one observed 
in Cas A I Laming and Hwanel lioOrO and in other 



remnant. However, given the lack of other emis- 
sion line in the spectrum and the overall weakness 
of the thermal component, a full comparison with 
the abundance patterns observed in ejecta of other 
young SNRs is not possible, so the interpretation 
of North Spur as a Si-rich ejecta knot needs more 
data to be confirmed. 

Assuming a reasonable line of sight of the 
emitting plasma in the North Spur (about 1 pc), 
a density of the knot of 3.8(1.6 — RS| v A ^™~ 3 



6) x d~ 



cm 



were derived (where da is the distan ce in units of 
5 kpc. lBecker and Szvmkowiaklll98lh . and a total 
mass of the thermal plasma of 0.7(0.3 — 1) X d§ 5 ^ 2 
M.q . Both density and mass are consistent with an 
ejecta of a Type II SN. However, the measured tem- 
perature (between 0.12 and 0.2 keV, considering 
both the EPIC and ACIS measurements) is lower 
then typical ejec ta knot temperature meas ured in 
young SNRs (e.g. lLaming and Hwanj20"03l for Cas 
A). 

This discrepancy can be eliminated if we 
assume that the ejecta we are observing have not 
been heated by the passage of the reverse shock, 
but by the supersonic shock of the expanding pul- 
sar wind nebula. The expansion of a PWN in the 
freely expanding eje cta of its SNR has been studied 
by several authors ^Reynolds and Cheyalierlll984l 
Ivan der Swaluw et alJl200ll IChevalierll2004 - The 
X-ray emission associated to the PWN supersonic 
shock into the ejecta have been obse rved up to now 
only in 3C58 l|Bocchino et alJl200lh . while for the 
Crab and 0540-69 has been observed in the optica l 
l|Sankrit and Hesterl Il997l iKirshner etafl Il989h . 
There are no other established observati onal signa- 
ture o f PWN expanding int o eiecta (seelChevalierl 
l2004l for a review). In 3C58. lBocchino et al.l(|200lT 
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measured an X-ray temperature of 0.2 keV, which 
is remarkably similar to the one found here in 
the North Spur of G21.5-0.9. However, the mass 
rep orted in the case of 3C58 was only 0.1 M Q 
and lChevalierl l|2004h noted that these was higher 
then the predicted swept-up mass of ejecta for 
the evolutionary stage of that remnant. If we as- 
sume a similar evolutionary stage for G21.5-0.9, 
our mass estimate in the North Spur is even more 
problematic for this scenario. 

We note that the interpretation of the North 
Spur in term of ejecta overrun by the PWN would 
imply a radius of the X-ray plerion much larger 
then the radio plerion (a bout 120 arcsec vs . 50 arc- 
sec of the radio plerion. iFurst et al1ll988h . a fea- 
ture which is not expected and indeed not observed 
in any other PWN. Even if such a large radius of 
the X-ray nebula compared to the radius of the 
radio nebula was suggested elsewher e in literature 
for G21. 5-0.9 l|Safi-Harb et alJl200lh . it would be 
hard to explain it in terms of conventional PWN 
models. 

Finally, we note that the differences of best- 
fit values between Chandra and XMM in Table 1 
do not affect the conclusion of the paper about the 
presence of an additional thermal component in the 
spectrum. One of the reason of the differences may 
be the different extraction radii used for Chandra 
and XMM (this fully explains the differences in the 
fluxes, for instance). Further discussions about the 
nature of the differences are outside the scope of 
this work. 
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